function fs = SGLD_features(image,d,sita)
%generate features based on SGLD
%param d: the ditance value of SGLD
%param sita: the angle value of SGLD
%generate four features:Energy,Entropy,Intensity,Inertia
sgldm = SGLD(image,d,sita);
%sgldm = sgldm/2000;
sgldm = sgldm*100/sum(sum(sgldm));
fs(1) = sum(sum((sgldm.*sgldm)));
fs(2) = sum(sum((sgldm.*log(sgldm+eps))));
[m,n] = size(sgldm);
inertia = 0;
intensity = 0;
for i=1:m
    for j=1:n
        intensity = intensity + sgldm(i,j)*i*j;
        inertia =  inertia + sgldm(i,j)*(i-j)*(i-j);
    end
end
fs(3) = intensity;
fs(4) = inertia;

